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We study the transport properties of an atomic-scale con- 
tact in the ballistic regime. The results for the conductance 
and related transmission eigenvalues show how the properties 
of the ideal semi-infinite leads (i.e. measuring device) as well 
as the coupling between the leads and the conductor influence 
the transport in a two-probe geometry. We observe the evo- 
lution from conductance quantization to resonant tunneling 
conductance peaks upon changing the hopping parameter in 
the disorder-free tight-binding Hamiltonian which describes 
the leads and the coupling to the sample. 

PACS numbers: 73.23.Ad 



Mesoscopic physicso has changed our understanding of 
transport in condensed matter systems. ilhe discovery 
of new effects, such as wak localization!] or universal 
conductance fluctuations^ has been accompanied by re- 
thinking of the established ideas in a new light. One of 
the most spectacular discoveries of mesoscopics is con- 
ductance quantization (CQ)Btl in a short and narrow 
constriction connecting two high-mobility (ballistic) two- 
dimensional electron gases. The conductance of these 
quantum point contacts as a function of the constriction 
width W ~ Xf has steps of magnitude 2e 2 /h. New ex- 
perimental techniques have allowed observation of simi- 
lar phenomenaQ in metallic point contacts of atomic size. 
The Landauer formula!] for the two-probe conductance 
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has provided an explanation of the stepwise conductance 
in terms of the number N < M of transverse propagat- 
ing states ( "channels" ) at the Fermi energy Ep which are 
populated in the constriction. Here t is the transmission 
matrix, T n transmission eigenvalues and Gq — 2e 2 /h is 
the conductance quantum. In the ballistic case (tt^)ij is 
5ij , or equivalently T n is 1 . 



• studies have explored 
-■They include geome- 
l1 ' temperature effects, 



Sij , or equivalently T n is 1 . Furthc 
CQ|-,uiidcr a range of conditions., 
try,EH13 scattering on impurities, 
and magnetic held. 

In this paper we study the influence of the attached 
leads on ballistic transport (£ > L, I being elastic mean 
free path, L being the system size) in a nanocrystal. We 
assume that in the two-probe theory an electron leav- 
ing the sample does not reenter the sample in a phase- 
coherent way. This means that at zero temperature phase 
coherence length is equal to the length of the sample 
L. In the jargon of quantum measurement theory, the 
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FIG. 1. A two dimensional version of our actual 3D model 
of a two-probe measuring geometry. Each site hosts a single 
s-orbital which hops to six (or fewer for surface atoms) near- 
est neighbors. The hopping matrix element is t (within the 
sample), £l (within the leads), and tc (coupling of the sam- 
ple to the leads). The leads are semi-infinite and connected 
at ±oo to reservoirs with potential difference [Xh — fJ-R = eV. 



leads act as a "macroscopic measurement apparatus". 
Our concern with the influence of the leads on conduc- 
tance is therefore also a concern of quantum measure- 
ment theory. Recently, the effects of a lead-sample con- 
tact on quantum transport in molecular devices have re- 
ceived increaseeLattention in the developing field of "na- 
noelectronics" £3 Also, the simplest lattice model and re- 
lated real-space Green function technique are chosen here 
in order to address some practical issues which appear in 
the frequent use of these methodsa to study transport 
in disordered samples. We emphasize that the relevant 
formulas for transport coefficients contain three different 
energy scales (corresponding to the lead, the sample, and 
the lead-sample contact), as discussed below. 

In order to isolate only these effects we pick the strip 
geometry in the two-probe measuring setup shown on 
Fig. 0. The nanocrystal ( "sample" ) is placed between two 
ideal (disorder-free) semi-infinite "leads" which are con- 
nected to macroscopic reservoirs. The electrochemical 
potential difference eV = hl — is measured between 
the reservoirs. The leads have the same cross section as 
the sample. This eliminates scattering induced by the 
wide to narrow geometryli3 of the sample-lead interface. 
The whole system is described by a clean tight-binding 
Hamiltonian (TBH) with nearest neighbor hopping pa- 
rameters i mn 
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FIG. 2. Conductance G of an atomic-scale ballistic con- 
tact 3x3x3 for the following values of lead and coupling 
parameters: (a) tc = 1,£l = 1, (b) tc = 1.5, £l = 1 (c) 
tc = 3, th = 1, and (d) tc = 0.1, tc = 1. In the case (d) 
the conductance peaks are connected by smooth curves of 
G < 0.004e 2 /fr. 

H=J2 *mn|m)(n|, (2) 

(m,n) 

where |m) is the orbital ^(r — m) on the site m. The 
"sample" is the central section with N x x N y x N z sites. 
The "sample" is perfectly ordered with t mn = t. The 
leads are the same except t mn = t^. Finally, the hopping 
parameter (coupling) between the sample and the lead 
is i mn = tc- We use hard wall boundary conditions in 
the y and z directions. The different hopping parameters 
introduced here have to be used to get the conductance 
at Fermi energies, throughout the whole band extended 
by the disorder ,£3 i.e. £l > t. Thus, one has to be aware 
of the conductances we calculate in our analysis when 
engaging in such studies. 

Our toy model shows exact conductance steps in multi- 
ples of Gq when tc = ti, = t. This is. a consequence of in- 
finitely smooth ("ideally adiabatic"EI) sample-lead geom- 
etry. Then we study the evolution of quantized conduc- 
tance into resonant tunneling conductance while chang- 
ing the parameter t^ of the leads as well as the coupling 
between the leads and the conductor tc- An example of 
this evolution is given on Fig. |^. The equivalent evolution 
of the transmission eigenvalues T n of channels is shown 
on Fig. |j| A similar evolution has been studied recently 
in one-atom point contacts.li-3 

The non-zero resistance is a purely geometrical ef- 
fccttL3 caused by reflection when the large number of 
channels in the macroscopic reservoirs matches the small 
number of channels in the lead. The sequence of steps 
(1, 3, 6, 5, 7, 5, 6, 3, 1 multiples of Gq as the Fermi energy 
Ep is varied) is explained as follows. The eigenstates 
in the leads, which comprise the scattering basis, have 
the form ip^ oc sin{k y m y ) sa\{k z m z )e lkxmx at atom m, 
with energy E = 2t\,[cos(k x a) + cos(k y a) + cos(fc 2 a)], 
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FIG. 3. Transmission eigenvalues of the atomic-scale bal- 
listic contact 3x3x3. The parameters £l and tc are the 
same as in Fig. |^. All channels = (k y (i), k z (j)) whose 

subbands are identical have the same T n . This gives the de- 
generacy of T n : three (upper panel), two (middle panel), and 
one (bottom panel). In the middle panel the lower two sub- 
bands have an energy interval of overlap with the upper two 
subbands. 

where a is the lattice constant. The discrete values 
k y (i) = in/(N y + l)a and k z (j) = j^/{N z + l)a de- 
fine subbands or "channels" labeled by (k y ,k z ) = 
where i runs from 1 to N y and j runs from 1 to N z . The 
channel (k v , k z ) is open if Ep lies between the bottom of 
the subband, 2tj\— 1 + cos(k y a) + cos(k z a)], and the top 
of the subband, 2*l[1 + cos(k y a) + cos(fc z a)]. Because of 
the degeneracy of different transverse modes in 3D, sev- 
eral channels (k y ,k z ) open or close at the same energy. 
Each channel contributes one conductance quantum Gq . 
This is shown on Fig. |^ for a sample with 3x3 cross sec- 
tion where the number of transverse propagating modes 
is M = 9. In the adiabatic geometry, channels do not 
mix, i.e. the transmission matrix is diagonal in the basis 
of channels defined by the leads. 

We compute the conductance using the expression ob- 
tained in the framework of Keldysh technique by treating 
the coupling between the central region and the lead as 
a perturbation.Ej This provides the following, Landauer- 
type, formula for the conductance in the non-interacting 
system 

G=^Tr(f L G r 1N J R G a Nxl ) =^Tr(ttt), (3) 

t = JF L G r 1N ^ R . (4) 

Here G r lN , G"^ 1 are matrices whose elements are the 
Green functions connecting the layer 1 and N x of the 
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sample. Thus only the block N y x N z of the whole matrix 
G(n, m) is needed to compute the conductance. The 
positive operator f L = i(t r L - ££) = -2ImE L > 
is the counterpart of the spectral function A = i(G r — 
G a ) for the self-energy Xj, introduced by the left lead. 
It "measures" the coupling of the open sample to the 
left lead (T^ is equivalent for the right lead). The Green 
operator is defined as the inverse of (E — H) including the 
relevant boundary conditions. Instead of inverting the 
infinite matrix we invert only [E — Hg) defined on tha 
Hilbert space spanned by orbitals |m) inside the sampled 



G' 



(E-H s - £ r y 



(5) 



where Hs is TBH for the sample only. This is achieved by 
using the retarded self energy Y7 — Y7 L + fip introduced 
by the left (L) and the right (R) lead. In site represen- 
tation Green operator G r ' a is a Green function matrix 
G r ' a (n, m) = (n|G"~' a |m). Equation (0) does not need 
the small imaginary part i0 + necessary to specify the 
boundary conditions for the retarded or advanced Green 
operator G r,a because the lead self-energy (X a = [S r ]T) 
adds a well defined imaginary part to E — H$- This 
imaginary part is related to the average time an electron 
spends inside the sample before escaping into the leads. 
The self-energy terms have non-zero matrix elements only 
on the edge lawrs of the sample adjacent to the leads. 
They are givenEJ in terms of the Green function on the 
lead edge layer and the coupling parameter tc 
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xt, r (k y , k z ) sm(k y m y ) sin(fc z m z ), (6) 

where (n, m) is the pair of sites on the surfaces inside the 
sample which are adjacent to the leads (L or R). The 
self-energy Y7(k y , k z ) in the channel (k z , k y ) is given by 
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for \E%\ < 2ti,. We use the shorthand notation E% = 
E — e(k y , k z ), where e(k y , k z ) — 2iL[cos(fc y a) + cos(fc z a)] 
is the energy of quantized transverse levels in the lead. 
In the opposite case \Ee\ > 2t^ we have 
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In order to study the conductance as a function of two 
parameters th and tc we change either one of them while 
holding the other fixed (at the unit of energy specified 
by t), or both at the same time. The first case is shown 
on Fig. 2 and Fig. |] (upper panel), while the second one 
on Fig. I (lower panel). The conductance is depressed in 
all cases since these configurations of hopping parameters 
tmn effectively act as a barriers. There is a reflection at 
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FIG. 4. Conductance G of the atomic-scale ballistic con- 
ductor 3x3x3 for the following values of lead and coupling 
parameters: Upper panel — (a) tc = 1, £l = 1, (b) tc = 1, 
tL ~ 1.5, and (c) tc = 1, th = 3; Lower panel — (a) tc = 1, 
th = 1, (b) tc — 1.5, tL = 1.5, and (c) tc = 3, th = 3. 



the sample- lead interface due to the mismatch of the sub- 
bands in the lead and in the sample when t-^ differs from 
t. This demonstrates that adiabaticity is not necessary 
condition for CQ. In the general case, each set of channels 
which have the same energy subband is characterized by 
its own transmission function T n (Ep). When the cou- 
pling tc = 0.1 is small a double-barrier structure is ob- 
tained which has a resonant tunneling conductance. The 
electron tunnels from one lead to the other via discrete 
eigenstates. The transmission function is composed of 
peaks centered at E r — 2i[cos(fc a: a)+cos(fc !/ a)+cos(/c z a)], 
where k x = kTr/(N x + l)a is now quantized inside the 
sample, i.e. k runs from 1 to N x . The magnitude and 
width of peaks is defined by the rate at which an electron 
placed between barriers leaks out into the lead. These 
rates are related to the level widths generated through 
the coupling to the leads. In our model they are en- 
ergy (i.e. mode) dependent. For example at Ep = 
seven transmission eigenvalues are non-zero (in accor- 
dance with open channels on Fig. |^) and exactly at 
Ep = three of them have T = 1 and four T = 0.5. 
Upon decreasing tc further all conductance peaks, ex- 
cept the one at Ep = 0, become negligible. Singular 
behavior of G(E&) at subband edges of the leads was 
observed before JIj 

It is worth mentioning that the same results are ob- 
tained ljtsing a non-standard version of Kubo- Greenwood 
formulaliU for the volume averaged conductance 
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G = — Tr hvjm G hv x lm G , (9a) 

h L x \ J 

ImG= (G r -G a ), (9b) 

where v x is the x component of the velocity operator. 
This was originally derived for an infinite system with- 
out any notion of leads and reservoirs. The crucial non- 
standard aspect is use of the Green function (|J) in for- 
mula (||). This takes into account, through lead self- 
energy (|6|), the boundary conditions at the reservoirs. 
The reservoirs are necessary in both Landauer and Kubo 
formulations of linear transport for open finite systems. 
They provide thermalization and thereby steady state pi 
the transport in the central region. Semi-infinite leadsES 
are a convenient method to model the macroscopic reser- 
voirs. When employing the Kubo formula ([)]) one can 
use current conservation and compute the trace only on 
two adjacent layers inside the sample. To get the correct, 
results in this scheme L x in Eq. (^) should be replacedES 
by a lattice constant a. 

In the quantum transport theory of disordered systems 
the influence of the leads on ±he conductance of the sam- 
ple is understood as follows O An isolated sample has a 
discrete energy spectrum. Attaching leads necessary for 
transport measurements will broaden energy levels. If the 
level width T due to the coupling to leads is larger than 
the Thouless energy E Th = h/r D ~ hD/L 2 , (D = v F £/3 
being the diffusion constant) the level discreteness is 
unimportant for transport. For our case of ballistic con- 
duction, i?Th is replaced by the inverse time of flight 
Hvf/L. In the disordered sample where T -&rh, vary- 
ing the strength of the coupling to the leads will not 
change the transport coefficients. In other words, the in- 
trinsic resistance of the sample is much larger than the 
resistance of the lead-sample contact. E£l In the opposite 
case, discreteness of levels becomes important and the 
strength of the coupling defines the conductance. This 
is the realm of quantum dotstil where weak enough cou- 
pling can make the charging energy e 2 /2C of a single 
electron important as well. Changing the properties of 
the dot-lead contact affects the conductance, i.e. the re- 
sult of measurement depends on the measuring process. 
The decay width T = ft/Vd W eii of the electron emission 
into one of the leads is determined by transmission prob- 
abilitieS|-of channels through the contact and mean level 
spacing.Ej This means that mean dwell time Tdwcii in- 
side our sample depends on both tc and t^. Changing 
the hopping parameters will make Tdweii greater than the 
time of flight r/ = L/vp. Thus we find that ballistic con- 
ductance sensitively depends on the parameters of the 
dephasing environment (i.e. the leads). 

In conclusion, we have studied the transport properties 
of a ballistic nanocrystal placed between two semi-infinite 
leads in the simplest strip geometry. We observe extreme 
sensitivity of the conductance to changes in the hopping 
parameter in the leads as well as the coupling between 
the leads and the sample. As can be easily anticipated, 



the conductance evolves from perfect quantization (as a 
result of an ideal adiabatic geometry) to resonant tun- 
neling. Nevertheless, it is quite amusing that vastly dif- 
ferent G(Ep) are obtained between these two limits (e.g. 
Fig. |^). The results are of relevance for the analogous 
theoretical studies in disordered conductors as well as in 
the experiments using clean metal junctions with differ- 
ent effective electron mass throughout the circuit. 
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